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Abstract. 

In a recent article, Alii and Sewell formulated a new version of the Dicke-Hepp-Lieb 
laser model in terms of quantum dynamical semigroups, and thereby extended the macro- 
scopic picture of the model. In the present article, we complement that picture with a 
corresponding microscopic one, which carries the following new results, (a) The local mi- 
croscopic dynamics of the model is piloted by the classical, macroscopic field, generated 
by the collective action of its components; (b) the global state of the system carries no 
correlations between its constituent atoms after transient effects have died out; and (c) in 
the latter situation, the state of the system at any time t maximises its entropy density, 
subject to the constraints imposed by the instantaneous values of its macroscopic variables. 
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I. Introduction 



The classic work by Hepp and Lieb^ 1 -* (HL) on the theory of the Dicke laser model 1 - 2 ) 
represents a landmark in the constructive quantum theory of phase transitions, far from 
equilibrium, in open quantum systems. Most notably, it demonstrated how the model un- 
dergoes a transition from normal to coherent radiation when the optical pumping strength 
reaches a critical value. This substantiated ideas that had been proposed by Graham and 
Haken^ 3 ) at a more heuristic level. 

In a recent article, which we shall refer to as (I), Alii and Sewell^ 4 ) constructed a 
new version^ 5 ) of the Dicke model, in which the dynamics of the matter and radiation is 
governed by a quantum dynamical semigroup, which incorporates the dissipative action of 
an array of pumps and sinks: by contrast, the HL model represented a conservative system, 
comprising the matter, radiation and these reservoirs. The principle new results of (I) were 

(a) a generalisation of the HL form of the dynamics of the matter-radiation composite, 

(b) a demonstration that the model exhibits optically chaotic phases, in addition to the 
normal and coherent ones, and (c) a generalisation of the theory of quantum dynamical 
semigroups'- 6 ) to ones with unbounded generators. 

The theories of both HL and (I) were centred on the macroscopic equations of motion, 
derived from the underlying quantum dynamics of their respective models in a limit where 
their sizes become infinite. In fact, they did not provide any treatment of the micrody- 
namics in this limit. The object of the present article is to complement the macroscopic 
picture of (I) with a corresponding microscopic one. 

In seeking a microscopic-cum-macroscopic picture of the dynamics of the model, we 
note that this has already been achieved for certain conservative systems with long range 
forces. For such systems, Morchio and Strocchi^ 7 ) have shown that the evolutes, b t , of 
the local microscopic observables, 6, depend not only on b and the time t, but also on 
certain macroscopic observables, i.e. 'observables at infinity': this represents a picture 
completely different from that of the Haag- Kastler scheme^ 8 -*, in which the dynamics of 
the micro- observables is autonomous. Furthermore, the explicit structure of the equations 
of motion for the microscopic and macroscopic observables, bt and B t , has been obtained 
by Bagarello and Morchio^ (BM) for a class of conservative, mean field theoretic models 
in the following form. 

§ f = *w (1.D 

Thus, the microscopic dynamics is piloted by the macroscopic observables, which, for their 
part, evolve according to a self-contained law. 

The question now arises as to whether a similar state of affairs prevails in the case of 
open dissipative systems, such as the laser model; and the first thing to say about that is 
that the methods of BM, which led to equn. (1.1), are not applicable there. The reason is 
simply that those methods depended heavily on the fact that, in a conservative system, the 
dynamics is automorphic and so preserves the algebraic structure of the observables. In an 
open dissipative system, on the other hand, the dynamics is governed by a semigroup^ 4 -' of 
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transformations of the observables that does not preserve that structure, i.e., in general, 
(ab) t ^a t b t . 

Thus, any attempt to derive equations of the form (1.1) for the laser model must be 
based on methods that are applicable to dissipative systems. Here, the essential problem is 
to pass to a limit where the size of the system becomes infinite, since the dynamics of the 
finite version of the model has already been formulated in (I). In fact, we surmount this 
problem by means of a compactness argument, that lends itself to the application of the 
Arzela-Ascoli theorem. As a result, we do indeed obtain a microscopic-cum-macroscopic 
dynamics of the form (1.1), though now, of course, the functions / and F incorporate the 
dissipative action of the pumps and sinks on the system. 

The principal new results that we obtain in this way are of a very simple character, 
and may be summarised as follows. 

(1) The microscopic observables of this particular model are those of the matter only, 
since the radiation field reduces to a classical macroscopic one. 

(2) The dynamics of the local microscopic observables is piloted by that classical 
radiation field. 

(3) The microdynamics contains certain transients, whose lifetimes are of the order of 
that of an atomic excited state. This dynamics simplifies considerably after these transients 
have died out, and leads to the decorrelation of the observables of different atoms. Further- 
more, the resultant global microstate of the system is stationary and space-translationally 
invariant in the normal phase, and periodic with respect to both space and time in the 
coherent one. Evidently, it is very complicated in the chaotic phase. 

(4) Following the decay of the transients, the instantaneous state of the model is 
determined by a non-equilibrium variational principle. Specifically, this state maximises 
the global entropy density, subject to the constraints imposed by the prevailing values 
of the macroscopic observables. As far as we know, this is the first example for which a 
principle of this kind has been proved for the non-equilibrium regime. 

We shall present our arguments as follows. In Sections II and III, we shall provide a 
brief summary of (I), re-expressed in a form suitable for our present purposes. Specifically, 
Section II will be devoted to a formulation of the model and Section III to a resume of 
its macroscopic properties. In Section IV, we shall derive the above results (1) and (2) 
by passing to the infinite size limit of the quantum dynamics of the finite version of the 
model. We shall then derive the results (3) and (4) in Sections V and VI, respectively. We 
shall conclude, in Section VII, with a brief discussion of our main results and their possible 
ramifications. 

II. The Model. 

Our model, as formulated in (I), is a dissipative quantum system, E^), consisting of 
a chain of 2N + 1 identical two-level atoms interacting with an n— mode radiation field. 
We label it by N only, since we shall be concerned with its properties in a limit where 
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N tends to infinity, while n remains fixed and finite. We formulate as a quantum 

dynamical system (B^ N \ T^ N \ «S^), where B^ is a algebra of observables, T^ N ^ is a 
one-parameter semigroup of completely positive (CP) transformations of B^ N \ and is 
a family of normal states on B^ N \ that is stable under the dual of T^ N \ Most importantly, 
T( N ) incorporates the dissipative action of pumps and sinks, which do not appear explicitly 
in the model, on the matter-radiation system. We build the model from its constituent 
parts as follows. 

The Single Atom. This is assumed to be a two-state atom or spin, T, at . Its algebra 
of observables, A a t, is that of the two-by-two matrices, and is thus the linear span of Pauli 
matrices (a x , a y , a z ) and the identity, /. Its algebraic structure is provided by the relations 

a l = a l = \ = I', a ^ a v = icr z, etc. (2.1) 
We define the spin raising and lowering operators 

o± = ^(a x ±ia y ). (2.2) 

We assume that the atom is coupled to a pump and a sink, and that accordingly (cf. (I)) 
its dynamics is given by a one-parameter semigroup {T at (t) |tGR+} of completely positive, 
identity preserving contractions of A a t, whose generator, L aU is of the following form. 

L a t<J± = -(7i=Be)cr±; L at a z = -7 2 (a 2 - rjl), (2.3) 

where e(> 0) is the energy difference between the ground and excited states of an atom, 
and the 7's and r\ are constants whose values are determined by the atomic coupling to 
the energy source and sink, and are subject to the restrictions that 

< 72<27i; -1<77<1. (2.4) 

In particular, r\ is positive or negative according to whether the coupling of the atom to 
the pump or the sink is the stronger. In the former case, these couplings drive the atom 
to a terminal mixed state with inverted population, i.e. with higher occupation probability 
for the excited state than for the ground state. 

The Matter. This consists of 2N + 1 non-interacting copies of E at , located at the 
sites r = — N, . ., iV of the one-dimensional lattice Z. Thus, at each site, r, there is a copy, 
S r , of E a t, whose algebra of observables, A r , and dynamical semigroup, T r , are isomorphic 
with A a t and T at , respectively. We denote by a r>u the copy of a u at r, for u = x,y, z, ±. 

We define the algebra of observables, A^ N \ and the dynamical semigroup, T^2t, of the 
matter to be ®^ = _ N A r and ®^ = _ N T r , respectively. Thus, ^4^) is the algebra of linear 
transformations of C 47V+2 . We identify elements A r of A r with those of A^ N > given by 
their tensor products with the identity operators attached to the remaining sites. Under 
this identification, the commutant, A' r , of A r is the tensor product ® s ^ r A s . 



4 



It follows from these specifications that the generator, L { m J t , of T^J. is given by the 
formula 

Lmat = ^r, (2.5) 

where 

L r 0y 5± = -(7i=Be)oY,±; i r ^, 2 = -l2{a r ,z - ^-f); 

and L r (A r A' r ) = (L r A r )A' r VA r eA r , A' r eA' r (2.6) 

The Radiation. We assume that the radiation field consists of n(< oo) modes, 
represented by creation and destruction operators {a*, ai\l = 0, . ., n — 1} in a Fock-Hilbert 
space Tirade as defined by the standard specifications that (a) these operators satisfy the 
CCR, 

[ai, a^} = Si m I; [a h a m ] = 0, (2.7) 

and (b) H ra d contains a (vacuum) vector $, that is annihilated by each of the a's and is 
cyclic w.r.t. the algebra of polynomials in the a*'s. 

We define the Weyl map z = (z , . ., z n -i)^W{z), of C n into C(H ra d) by the formula 

En— 1 
Z;a; (2.8) 

We then define 7?. to be the *— algebra of polynomials in the a's, a*'s and the Weyl operators 
W{z), with z running through C n . Thus, in view of the CCR (2.7), this algebra is the linear 
span of terms of the form W(z)a\ . .a\ , where a* denotes a or a*. Equivalently, it is the 
linear span of the derivatives, of all orders, of the operators W(toZo, . .,t n -iz n -i) w.r.t. 
the real variables (to, • -,t n -i). 

We assume that the radiation dynamics is given by the canonical extension to 1Z of 
Vanheuverszwijn semigroup^ 10 -*, T ra ^, of quasi-free CP transformations of the Weyl algebra 
of linear combinations of {W (z)\z<EC n } . The formal generator of this semigroup is 

En— 1 
l=Q {m[a~tai, .] + 2« I a£(.)aj - Ki{a*a h .}), (2.9) 

where {., .} denotes anticommutator, and the frequencies, u>i, and the damping constants, 
ki, are positive. 

The Composite System. This is the coupled system, Y,( N \ comprising the matter 
and radiation. We assume that its algebra of observables, B^ N \ is the tensor product 
A {N) ®11. Thus, BW, like K, is an algebra of both bounded and unbounded operators in 
the Hilbert space H w := C 4N+2 ®H rad . We shall identify elements A, R, of A w , TZ, 
with A®I ra d and I mat <S>R, respectively. 

We assume that the matter-radiation coupling is dipolar and is given by the interaction 
Hamiltonian 

H™ = i{2N + l)-V2^_^^ = - 1 Az((Tr) _ arexp( _2 7rz / r /n) - h.c), (2.10) 
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where the A's are real-valued, N— independent coupling constants. Correspondingly, we 
define the radiation field, (p^, of the model by the stipulation that its value at the site r 
is the coefficient of a T}+ in this formula. Thus, 

n— 1 

0(^0 = -i(2N + l)~ 1/2 J2i- X i a i e M^ilr/n). (2.11) 



We now need to define the state space, S^ N \ and the dynamical semigroup, T^ N \ in a 
way that takes account of the unboundedness both of H^t and of some of the elements of 
B^. To this end, we start by defining "Z^o to be the space of density matrices in TL^ N \ 
with the trace norm topology, and T>[ to be the subset of its elements, p^ N \ for which 
Tr(B* p( N ) B) is finite for all B in B^ N ^ . We define Si to be the set of positive, normalised, 
linear functionals, V 7 , on B ■, that are in the one-to-one correspondence with the 
v[ N — class density matrices, given by the standard formula ^ N \B) = Tr(p( N ^B). We 
then adopt the formulation of and provided by (I), in which is constructed 
as the modification of T^ t ®T rad due to the interaction H^ t \ and is a 'large' subset 

of S{ N ' that is stable w.r.t. this dynamical semigroup. Specifically, the construction yields 
the following results^ . 

(1) The density matrices corresponding to form a dense subset of Vq N \ in the 
topology corresponding to the trace norm. 

(2) is stable under the transformations ^ N ^iP^oT^ N \t) := ^ t (A ° . 
(3) 

^-4 N) (B) = 4 N) (L^B) VBeBW tGR+, (2.12) 

where 

LW = LW + L rad + i[HS$,.]. (2.13) 
This completes our specification of E^ = (6^,^,5^). 

The Sequence of Systems {E^}. Since we shall be concerned with the properties 
of the above model in a limit where N tends to infinity and n remains fixed, we need 
to consider the sequence of systems {E^^ATeN}. Specifically, we shall investigate the 
evolution of these systems from initial states, ip( N \ for which the mean photon number is 
not super-extensive, i.e., its ratio to N is bounded. Thus, we assume that, for some finite 
constant C, 

ip {N \a^ai) < CN, ViVeN, Ze[0, n - 1]. (2.14) 

This exclusion of super-extensivity was shown in (I) to persist, in a uniform way, at later 
times, in that 

V> t (JV) (a?aj) < DN, ViVeN, £eR+, le[0, n - 1], (2.15) 
where D is a finite constant. 
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III. Macroscopic Dynamics and Phase Structure 

We formulate the macroscopic description of the model, as in (I), in terms of the 
global intensive observables 

s[ N) = (2N + l)" 1 ^ N __ Ar ^r _exp(-27rz/r/n); 1 = 0,.. .,n-l (3.1) 

and 

p\ N) = (2N + l)- 1 ^ = _ iv a r , 2 exp(-27rz/r/n); I = 0, .. .,n - 1, (3.2) 
together with the operators 

a\ N) = (2N + l)" 1/2 a z ; I = 0, .. .,n - 1, (3.3) 

corresponding to a scaling of the number operators a^cii in units of 2N + 1. We note that 
the set of p( N > : s is the same as that of their adjoints, since 

p™* = P ( N) ; and P W* = pl N _) for I = 1, . ., n - 1. (3.4) 

The set M^ N ^ := {sj^\sj^*,pj , , } of macroscopic variables is a Lie algebra 
w.r.t. commutation, and, by equns. (2.10), (2.11) and (3.3), the radiation field and the 
interaction Hamiltonian are the algebraic functions of them given by 

= -iY^~^Xia < l N) exp(2Tvirl/n) (3.5) 

and 

i4i> = + D5:;;;A,( S r»ar»* - ft . c .). 0.6) 

In (I), the dynamics of Af( N ) was extracted from the microscopic equation of motion 
(2.12), in a limit where N-^oo and n remains fixed and finite, subject to the assumptions 
of Section 2. As one might anticipate from the fact that the commutators of the elements 
of tend normwise to zero in this limit, this macrodynamics is classical. 

To describe this dynamics, we simplify our notation by defining x^" 1 = (x\ N \ . ., 
to be the 5n self-adjoint operators, which, by equns. (3.1)-(3.4), comprise the Hermitian 
and anti-Hermitian parts of the elements of M^ N \ Correspondingly, in anticipation of a 
classical limit for the dynamics of these observables, as iV^oo, we introduce a phase space, 
X = R 5n and, for each point x = (xi, . .,X5 n )&X, we define si,pi\l = 0, . .,n — 1} to 
be the complex numbers related to x in precisely the same way that the elements of M^ N ^ 
are to x^ N \ We define Cq(X) to be the space of bounded continuous functions on X, that 
tend to zero at infinity. 

We formulate the macroscopic dynamics of in terms of the time-dependent quan- 
tum characteristic function fi[ N ^ on X, defined by the formula 

4 N \v) = 4 N) {eMw-x {N) )) VyeX,teR + . (3.7) 
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Further (cf. (I)), we relate the limiting form of this function, as iV^oo, to the following 
classical dynamical system. 

Definition 3.1. We define /C to be the flow in X, whose equation of motion, 

£ - FM, (3.8) 

is of the explicit form 

= -{iui + Ki)ai tt + Xisi >t (3.9a) 

E 1 n—l 
X m P\l-m],t Oi m t yo.tlO) 
m=0 1 J 



dt 
and 

= ~l2(pi,t - vho) - 2^ m=o A m (a mit s [i+m])t + a m ,tS[ m -i],t)> ( 3 - 9c ) 
where [l±m] — l±m (mod n). 

The following Propositions were proved in (I). 

Proposition 3.2. The equation of motion (3.8) (i.e. (3.9)) has a unique global 
solution, corresponding to a one- parameter semigroup of transformations x—>Ttx(=xt) of 
X, that maps the compacts into compacts. 

Proposition 3.3. (%) Under the above assumptions, /i[ N ^ tends pointwise, as iV^oo, 
to the characteristic function, \it, of a probability measure, mt, on X, i.e. 

limjv^ooAtt (y) = Ih{v) = / expi(x.y)dm t (x), (3.10) 

J x 

convergence being uniform on the X — compacts. 

(ii) The evolution of m t is induced by the dynamical semigroup, r, of K, according to 
the formula 

[ dm t f= [ dm Q for(t) V/gC (X), teR+. (3.11) 
J x J x 

(Hi) In the particular case of a pure phase, where mo is a Dirac measure, 8 X , with 
support at some point x, then m t = S Xt . 

Comment. This last Proposition signifies that, in the limit iV^oo, the macroscopic 
variables (a^ N \ s^ N \p^) reduce to classical ones (a, s,p). Correspondingly, the radiative 
field, <f>( N \ as given by equn. (3.5), reduces to a classical one, cf>, whose value at the site r 
and time t is 

En — 1 
;=Q A/oj t) /exp(27uW/n). (3.12) 
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Phase Structure. By Def. 3.1, Prop. 3.2 and Prop. 3.3(iii), the properties of the 
pure phase are governed by the classical dynamical system /C. A study of this system yields 
the following results (cf. (I)). 

(1) For rj less than a certain specified critical value, rji (> 0), the system has a unique 
stable stationary state corresponding to the fixed point, xo(gX), given by 

ai =0; s t = 0; Pl = rjS l0 V/g[0, n - 1]. (3.13) 

(2) As rj increases through a certain positive value, 771 , the fixed point xq becomes 
unstable, and, by Hopf bifurcation, gives way to a periodic orbit of the form 

a ht = a z (0) expH(z4 + 6^))S lk ; s M = sj 0) exp((-^ + 9^))S lk ; p M = Vl 6 l0 , (3.14) 

where the selection of the /c'th mode and the values of the constants a^°\ and v are 
determined by the parameters of the model, while the phase angle 9^ is indeterminate. 
Thus, in optical terms, there is a transition from normal to coherent radiation as rj passes 
through 771 . This entails a breakdown of the gauge symmetry, represented by the trans- 
formations in which the cr rj _,az, and correspondingly are all rephased by the same 
factor exp(z6>). 

(3) The model generically undergoes a further transition from coherent to chaotic 
radiation, as rj passes through a second critical value, rj 2 (> 771). Specifically, the following 
two scenarios, both of which involve gauge symmetry breakdown, are feasible. 

(a) There can be a bifurcation, at rj = 772, from a periodic orbit to a strange attractor, 
as specified in the Ruelle-Takens scheme^ 11 ). This scenario is actually realised in the single 
mode case^. 

(b) For large n, there could be a succession of bifurcations, corresponding to the 
activation of different modes, leading to a chaotic phase along the lines of Landau's theory 
of hydrodynamical turbulence^ 12 -*. 

IV. The Microscopic Dynamics in Pure Phases. 

Our objective now is to obtain the microscopic dynamics of the model, in the limit 
N^oo. We take the local microscopic observables to be those of the matter only, since, as 
noted in the Comment following Prop. 3.3, the radiation field reduces to a classical macro- 
scopic one in this limit, and therefore this field carries no local microscopic observables 
within the terms of this model^ 13 ). 

We formulate the C*— algebra of the observables of the matter, in the limit N-^oo, as 
the standard quasi-local one for the infinite chain of atoms of the given species, occupying 
the sites of the lattice Z. Thus, in order to construct this algebra, we define L to be 
the family of all finite point subsets of Z, and we assign to each AgL the C*— algebra 
Aa := ®reA-4r, with A r as defined in Section 2. For AcA', we identify elements A of Aa 
with A®Ia'\a (€-Aa'), thereby rendering Aa isotonic w.r.t. A. We then define A to be the 
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completion of the local algebra Al '■= Ua6L^ a w - r -t- the norm it inherits from the Aa's. 
Thus, the algebras A^ N ' of the material observables of the systems E^ are just the local 
subalgebras A[-n,n] of A. 

We denote by S(A) the set of all states on A, and we represent the space translation 
group Z by the automorphisms V of this algebra, defined by the formula 

V(m)a rjU = a r +m,u, Vr, mGZ, u = ±,z. (4.1) 

We note that it follows from our specifications of the local structure of A that the 
semigroup T^ t , defined in Section 2, is just the restriction to A^ of the global one- 
parameter semigroup T mat := ®T r of CP transformations of A. The generator of T mat is 
therefore 

Lmat = / L r , (4-2) 

where L r is given by equn. (2.6). The effects of the matter-field coupling on the properties 
of Y,( N \ in the limit N^oo, will be treated presently. 

Our aim now will be to obtain the microdynamics of the matter induced by that of 
the systems Y*( N ' in the limit N^oo, subject to the following conditions, characteristic of 
pure phases, on their initial states. 

(PI) The restrictions of the initial states, ip( N \ of the systems E^\ to the algebras 
—n,n] are the same as those of a certain primary element, u, of S(A), i.e., 

(P2) cu(s| Ar ' > ) and u{pi) converge to values si and pi, respectively, for each /g[0, n— 1], 
as N^oo. 

(P3) The mean and dispersion of for the state ifj( N \ converge to values ol\ and 

zero, respectively, for each Ze[0, n — 1], as N^oo. 

The following lemma serves to confirm that the assumptions (PI- 3), together with 
equn. (3.7), imply that the conditions of Prop. 3.3(iii) prevail. 

Lemma 4.1. Under the above assumptions, the macroscopic probability measure, nit, 
is simply S Xt , where x t = r(t)x and x is the initial phase point corresponding to (a,s,p) 
of conditions (P2,3). 

Proof. By equn. (3.7) and the Schwartz inequality, 

|^ (y) - expiiy.x) N|^ (iV) (eMiy-(x (N) ~ *)) - I) |<2 (sin 2 (^^))] 1/2 

<l^{y.{x^-x)f] 1/2 , 
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and, in view of (P2,3), this tends to zero as N^oo. Hence, by Prop. 3.3(i), Ho(y) = 
exp(iy.x), i.e., mo = S x ; and consequently, by Prop. 3.3(iii), m t = S Xt , as required. 

The following Proposition signifies that, in the limit iV^oo, the microscopic dynamics 
is precisely that of the matter, under the action of the classical macroscopic radiation field, 

cf>. 

Proposition 4.2. (i) Assuming the conditions (Pl-3), the microscopic dynamics of 
the model corresponds, in the limit iV^oo, to a two-parameter family, 

{T(s, t\<p)\0<s<t; T{s, u\<t>)%{u, t\(f>) = T{s, t\<f>)}, 
of CP contractions of A, in that 

lim Jv ^ooV' (JV) (T {N) (t)A) = oo (T(0, t\<j>)A) VAg^l, teK+. (4.3) 
Furthermore, the generator of T is 

d 

£(t\</>) := —T(s,t\(j)) ls=t = L mat + 5^ r€Z [0r,tffr,+ - h.c, .]. (4.4) 

Thus, 

T(s,£|0) = T a exp( / duC{u\<t>)), (4.5) 

J s 

where T a is the antichronological operator. We define 

uj t :=ojoT(0,t). (4.6) 

(ii) T factorises into a product of single site contributions according to the formula 

T(s,t\(j)) = ^ r&z T r (s,t\(j) r ), (4.7) 



where 



T r (s,t\(f> r ) = T a exp( I duC r {u\<i>)), (4.8) 

J s 



and 

C r {t\4>) = L r + [(p r ,to- r ,+ - h.c, .]. (4.9) 



Comment. This result is similar to that obtained by Bagarello and Morchio^ 9 '' for a 
class of conservative systems, in that the microscopic evolution is governed by a classical 
field, formed by the time-dependent macroscopic observables. 



The proof of the Proposition requires the following lemmas. 
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Lemma 4.3. Defining 

a\f := a\ N ^ ~ a l|t , (4.10) 

^fC^i?) = lim^(ajffiW*) = 0, (4.11) 
the convergence being uniform on the t— compacts. 

Lemma 4.4. Defining 



z(N) ._ (N) r419 x 



lim^ oo ^ (JV) (4f4f*) = 0, (4.13) 
the convergence being uniform on the t-compacts. 

Proof of Prop. 4.2, assuming Lemma 4.3. (i) By equns. (2.5), (2.6), (2.9), 
(2.10), (2.12) and (2.13), 



f 

Jo 



4 N \A)-4 N \A)= [ ds^ s N \L mat A)- 

Jo 



where 



6? z = ±\ l Y* exp(±27uW/n)[(7 r± ,.]. (4.15) 

We note now that Al is stable under 5^, and that, by equns. (2.14), (2.15) and (3.3), 
a^) is bounded, uniformly w.r.t. N and t. It therefore follows from the 
Schwartz inequality that, for each A^Al, the norm of the r.h.s. of equn. (4.14) is likewise 
uniformly bounded. Hence, by the Arzela-Ascoli theorem, ^>[ N \A) converges to a limit, 
uniformly over the t— compacts, as N^oo over some subsequence of N. Therefore, as A 
is norm-separable, it follows from the standard diagonalisation procedure that there is a 
state u' t on A, such that, over some subsequence of the integers, 

hniA^ooV^A) = u' t (A) \/AeA L ,teTL+. (4.16) 

Further, by equns. (2.7), (4.2), (4.4), (4.8) and (4.10), we may re-express equn. (4.14) in 
the form 

(4.17) 

By Schwartz's inequality and Lemma 4.3, the last integral in this equation tends to zero, 
uniformly on the compacts, as iV— >oo. Hence, by equn. (4.16), we see that equn. (4.17) 
reduces to the following form in this limit. 



U ' t (A) - u/ (A) = f dsJ s {C{s\<j>)A), 
Jo 
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and hence 

^u/ t (A)=u/ t (C(tm). (4.18) 
Consequently, since equn. (4.5) implies that 

it follows from equn. (4.18) that 

d 

—(j'(T(8,t\(j>)A) = V0<s<t, AeA L . 

OS 

Hence, uj' s (T(s, t\(f))Aj is independent of s, and consequently, by assumption (PI) and 
equn. (4.16), 

u' t {A) = u f o (T(0,t\ct>)A)=u>(T(0,t\ ( f>)A), 

i.e., by equn. (4.6), 

u' t =u t =uoT(0,t\(j)). (4.19) 

Equations (4.3) and (4.4) now follow, at least for convergence over a subsequence of inte- 
gers, from equns. (4.16), (4.18) and (4.19). The removal of the restriction to a subsequence 
is achieved by noting that, in view of the uniqueness of the evolution ui— >u> t , specified by 
equn. (4.6), the compactness argument underlying the derivation of equn. (4.19) may be 
re-employed to show that all subsequential limits of {ijj( N \T( N>) (t)A)\NeN} are equal to 
u t (A). 

(ii) The factorisation property is an immediate consequence of the fact that, by equns. 
(4.4) and (4.9), 

Proof of Lemma 4.3, assuming Lemma 4.4. Since, by equns. (2.7), (3.3) and 

(4.10), [oifl , Oifl ] = AT -1 , it suffices to establish that the l.h.s. of equn. (4.11) tends to 
zero, uniformly on the t— compacts, as N^oo. 

To this end, we define 

and then infer from this equation and equns. (2.12) and (4.10) that 

rl tW « 

aj l.t <(N)/,t(m\ . O 



^ = ^»((L<"> + M*^>). 



dt Tt vv dt 

It follows from this equation, together with equns. (2.13), (3.1)-(3.3), (4.10), (4.12) and 
(4.20) that 
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and consequently, 

j/f = j/^exp(-2M) + 2A,ite T dt , exp(-2K I (< " t'))^\a ( l N t , ) * sfi)). (4.22) 
Therefore, defining 

it follows from the Schwartz inequality and equn. (4.22) that 

4P<4" ) eM-^it)+ fdt'lMeM^tit-mK^J^] 1 / 2 . (4.24) 

Jo 

For any fixed T(> 0), we now define 

:= sup te[0 , T] 4P (4.25) 

and 

£><?? := sup tmT] K^\ (4.26) 

noting that these are both finite, in view of the equn. (2.15) and the uniform boundedness 
of «W. Hence, by equns. (4.24)-(4.26), 



Ju ) <J^ ) exp(-2^t) + 2|Aj| / ^ / exp(-2«/(t - t / ))[C , I ( ? ) D, ( ^ ) ] 1/2 V te[0,T}. (4.27) 

From this formula it follows easily that 

T (N) < AN) |A/| (at) n (tf)ii/2 



and hence that 

I A; 



^<^ + ^[^S?] 1/a , (4-28) 



which in turn implies that 



C<?»<i[M (BS ) ) i/ 2 + ( | («) +4 ^))V 2) ]^. (429) 

Further, by definition (4.20) and assumption (P3), jffl converges to zero as N^oo; and 

by Lemma 4.4, the same is true of . Hence, by equn. (4.29), cffl tends to zero as 
N^oo. By equns. (4.20) and (4.25), this is just the required result. 

Proof of Lemma 4.4. Let qj N ^ denote the Hermitian or anti-Hermitian part of s\ N ^; 
and, correspondingly, let q^ t denote the real or imaginary part of s^ t - Then, defining 

c\ N t \v) := ^ (A °(exp(^ (A °)) VveR, (4.30) 
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we see from equn. (3.7) that cf^ is just the restriction of the characteristic function n\ N ^ 
to one matter mode. Thus, by Prop. 3.3(iii) and Lemma 4.1, 

\\m. N ^ 00 c\ N t \v) = exp(ivq t ,i) , \/veR,teR+, (4-31) 

convergence being uniform on the compacts. On integrating this equation against the 
Fourier transform, f(v), of any element, /, of the Schwartz space V(R), we see that 

hm^ 0O ^ (Ar) (/( 9 ; 7V) )) = f( %t ) V/GT?(R),teR+, (4.32) 

convergence being uniform on the t— compacts for given /. We shall now use this formula 
to prove the required result, in the form of the equation 

hm N ^ t N \( q l N Y) = (Qi,t) m , for m = 1,2, (4.33) 

convergence being uniform on the compacts. 

To establish this formula, we first note that, by equns. (2.2) and (3.1), together with 
our definition of q$ N \ the norm of this operator is 1/2. Thus, if / is any X>(R)— class 
function with support disjoint from [—1/2,1/2], then /(<?} ) = 0, and hence, by equn. 
(4.32), f(qi,t) = 0. This implies that qij must lie in the interval [—1/2,1/2]. We now re- 
employ equn. (4.32), choosing / to be of the form f(q) = q m h(q), where m = 1 or 2, and 
h is an element of V(R), whose value is unity over the interval [—1/2,1/2]. This choice 
yields precisely the required result (4.33). 

V. Asymptotic form of the Dynamics and States of Pure Phases 

We shall now show that the microscopic dynamics of the model simplifies greatly when 
transient effects, which decay in times of the order of the lifetime of the excited states of 
its atoms, are discarded. 

We start by noting that, by Prop. 4.2, the microscopic dynamics of the model is 
completely determined by the action of the CP contractions, T r , on the single site algebras 
A r . Thus, defining 

<r r ,u(t) := T r (0, t)a r ,u, for u = +,-, x, y, z, (5.1) 
we see from equations (2.1), (4.5) and (4.7)-(4.9) that 

da r -(t) 



dt 

and 



+ (ie + 7i)cr r _(£) - i(j) r ^a r ^{t) = (5.2) 



<}a '' ={t) • -,.rr,, : (/) 2/(o„ / o„ (/) h.O --,///. (5.3) 



dt 
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Thus, the spin vector a r (t) = (a rjX (t), cr rj/ (£), cr r>z (t)) evolves according to a linear inho- 
mogeneous equation 

^- = b r (t)a r (t) + cl, (5.4) 

where b r (t) is a linear transformation of C 3 , whose explicit form is obtained from the 
equivalence of equn. (5.4) to the pair of equations (5.2) and (5.3), and 

c= (0,0,72??). (5.5) 
The solution of equn. (5.4) may be expressed in the form 

<7r(t) = 9r(t,0)a r (0) + / dt'g r (t,t')cl, (5.6) 

Jo 

where g r is the Green function for b r , determined by the formula 

^g r (t, t') = b r (t)g r (t, t') Vt>t'>0; g r (t, t) = 7 3 , (5.7) 

and I3 is the identity operator in C 3 . 

Lemma 5.1. For fixed, t', g(t,t') decays to zero, with increasing t, at least as fast as 
exp(— 7t), where 7 = 771171(71,72). 

Proof. Let £ r be an arbitrary element of R 3 , and, for fixed f (eR+) and any t > t', 
let £ r (t) = (&■,*(*), £r,y (*)>&■,*(*)) := 9r{t,t')£ r . Then, by equn. (5.7), £ r (t) satisfies the 
homogeneous linear equation obtained by replacing o~ r by £ r in equn. (5.4) and discarding 
the term cl on the r.h.s. Hence, defining £ rj ± = l/2(£ riX ±i£ rjy ), it follows from the 
equivalence between equn. (5.4) and the pair of equations (5.2) and (5.3) that the equation 
of motion for £ r (t) is also given by replacing a by £ in the latter equations and discarding 
the r.h.s. of (5.3). Thus, 

d 

fotr,-(t) + (71 - ie)f r - i(pr,t£r,z(t) = 

and 

jfrAt) +I2&,z{t) ~ 2i(0 r)t £r,-(*) " C.C.) = 0. 

Hence, 

^\tr(t)\ 2 =^(*\tr,-(t)\* + tr,*(t)*) = -8^1 |Cr,-(*) | 2 - ^^{tf = 

-2[7 1 £ riX (t) 2 + 7 1 ^(t) 2 + 7 2 £ r)Z (t) 2 ]<-2mzn( 7l ,7 2 )|£ r (t)| 2 . 
Consequently, defining 7 := 771171(71,72), 

|(|£ r (t)| 2 exp(2 7 t))<0, 
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which implies that |£ r (*)|> an d hence g r (t,t'), decay to zero at least as fast as exp(— 7^), 
with increasing t. 

Comment. This lemma signifies that the first term on the r.h.s. of equn. (5.6) is 
merely a transient, which decays in a time, 7 -1 , of the order of the lifetime of an atomic 
excited state. Consequently, by (5.5), the asymptotic form for <r r (t), obtained by discarding 
the transients, is 

a r (t) = O r (t)I, (5.8) 
where ^ 

9 r (t) = [ dt'g r (t,t')(0,0,<y 2V ). (5.9) 
Jo 

Hence, since the state u t is given by its action on the finite monomials in the components 
of the cr r 's, the following Proposition is a simple consequence of these last two equations, 
together with Prop. 4.1. 

Proposition 5.2. The asymptotic form of the time- dependent state u t carries no 
correlations between the observables at the different sites and is therefore of the form 

U t = ®r6ZW r ,i, (5.10) 

where the atomic state u r j is given by 

Mr,tM = 6 r (t). (5.11) 



Comment. It follows immediately from this Proposition that the asymptotic state 
of the model is completely determined by the form of 6 r . We shall denote the components 
of 6 r analogously with those of a r : thus 9 r = (# r)X , @r,y, 6r,z) an d 9 r ,± = \{dr,x^r,y)- The 
following Proposition provides an explicit formula for 9 r in terms of the asymptotic forms 
of the macroscopic variables s and p. 

Proposition 5.3. The function 6 r is given by the equations 

6> r ,_(t) = Y, l=0 8i, t exp(2mrl/n) (5.12) 

and 

r A 1 ) = 2^i =0 Pi,texp(2mrl/n). (5.13) 

Proof. By equns. (5.2), (5.3), (5.7) and (5.9), the equations of motion for 9 r are 
formally the same as those for oy, i.e., 

d6 (t) 

r ' df { ) + (ie + -yi)O r ,-(t)-i<l> r ,tO r , z (t) = (5.14) 

and 

+ " 2i(0 r>t r ,-(f) - c.c.) = l2V . (5.15) 



17 



Further, denning 



En — 1 
^ s; ) texp(27rir//n) (5.16) 



and 

Pr(*) = 2^ =0 Putex-V&nirl/n), (5.17) 
we see, by elementary Fourier analysis, that equns. (3.9b) and (3.9c) are equivalent to 

+ (ie + 7l)«r(*) " i0r,tPr(*) = (5.18) 

and 

civ (t) — 

+ 72P r (t) - 2i{(j> r /s r {t) - c.c.) = 72 77, (5.19) 

respectively. Hence, the equations of motion (5.14) and (5.15) for Q r _ and 9 rjZ are equiv- 
alent to those given by (5.18) and (5.19) for s r and p r . Consequently, in the asymptotic 
regime, where the initial conditions are irrelevant, 6> r ,_ = s r and 6> r , 2 = p r . By equations 
(5.16) and (5.17), this is the required result. 

The following Corollary is an immediate consequence of Props. 5.2 and 5.3. 

Corollary 5.4. The asymptotic state ui t is spatially periodic, with periodicity n, i.e., 

u t = w t oV(n), (5.20) 

where V is the group of spatial translational automorphisms of A defined by equn. (4-1)- 

Explicit Asymptotic Form of ut- Equations (5.10)-(5.13) serve to express the time- 
dependent microstate Ut in terms of the macroscopic variables. We shall now present its 
explicit form for both the normal and coherent phases. We have no means of doing the 
same thing for the chaotic phase, since we have no explicit solution for the macroscopic 
dynamics there. 

1. The Normal Phase. Here, by equns. (3.13) and (5.10)- (5.13), u t takes the 
stationary asymptotic value Wqo, defined by the equations 

Uoo = ®r£ZWoo,r (5.21a) 

and 

Woo,r(ffr,±) = 0; Woo,r(o-r,«) = V- (5.216) 

2. The Coherent Phase. In this phase, it follows from equns. (3.14) and (5.10)- 
(5.13) that the asymptotic form, u)l° h ', of uj t is given by the formula 

w? h = ®r € zw£?, (5.22a) 

where 

= *! 0) exp(z(^ -ut- *„)); <tK.) = V- (5-226) 
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Comment. By contrast with the thermal equilibrium states of pure phases, the stable 
states (jj1° h are time-dependent. The same is true for the pure chaotic phases, since the 
macrosopic variables s t and p t are time-dependent there too* 14 ). 



VI. A Non-Equilibrium Maximum Entropy Principle 

We shall now establish that, as a further consequence of the results of Section 5, ut 
maximises the global entropy density of the matter, as defined on the states with spatial 
periodicity n, subject to the condition that the limiting value, as A/— >oo, of the expectation 
values of the macroscopic variables s( N \p( N "> are s t ,Pt, respectively. 

For this purpose, we start by defining S n to be the set of spatially periodic states, tp, 
on A, with period n, and <S n ,x t to be the subset of these states for which 

lim A r^ 00 'i/'(s[ Ar:) ) = s/,* and \im N ^ 00 ifj(p < i N) ) = pi, t , V7e[0, n], £eR+. 

Equivalently, by equns. (3.1) and (3.2), S HjXt consists of the states, ip, of spatial periodicity 
n, which satisfy the conditions 

— i s^ n ~^ 

n 2^ r=0 ^( cr r,-)exp(-27rz/r/n) = s Z)t (6.1) 

and 

n /> ^V'(^r,z)exp(-27rz/r/?i) = pi, t , (6.2) 

'r=U 

for / = 0, . ., n — 1 and tGR+. 

We formulate the global entropy density functional, s on S n in the following standard 
way^ 15 ). For ipES n and AgL, we define p\ to be the density matrix in 7i(A) representing 
the restriction of ip to ^4(A), and S\(ip) to be the corresponding local entropy, i.e., 

S A (V0 = -kTr(ptlnpt), (6.3) 

k being Boltzmann's constant. We define the global entropy density functional, s, on S n 
by the formula 

s(V>) = lim ATZ » ( 6 - 4 ) 

the convergence being guaranteed by the subadditivity of entropy* 16 ) and the periodicity 
of ip. 

Proposition 6.1. The asymptotic time- dependent state ui t maximises the functional 
s, as restricted to S njXt . 

Proof. We first note that, by Cor. 5.4, uj t eS n7 and that, by equations (5.10)-(5.13), 
it satisfies the conditions (6.1) and (6.2). Hence, u t ^S njXt . 

Next, we observe that, if ip is an arbitrary element of S njXt , then it follows by ele- 
mentary Fourier analysis from equns. (5.12), (5.13), (6.1) and (6.2) that ip(a r ) = 9 r (t), 
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for re[0, n — 1], and consequently, by the periodicity of tp, for all rGZ. Therefore, by equn. 
(5.11), the restrictions of ip and u>t to each single site algebra A r are identical. In other 
words, all the elements of S njXt coincide on each of the single site algebras. 

Further, by the subadditivity of entropy, 

< 'r 6 A 1 ' 

while, by equns. (5.10) and (6.3), 

S A (u t ) = J2 r€A S {r}M- 

Consequently, as ip and uj t coincide on the single site algebras, 
and therefore, by equn. (6.4), 

s(ip)<s(u) t ) VipeS n}Xt , 

which proves the Proposition. 

VII. Concluding Remarks. 

We have provided a microscopic picture of the multi-mode laser model, that is com- 
plementary to the macroscopic one of (I). Our principal new results are the following. 

(1) The microscopic dynamics is governed by the classical, macroscopic field, 0, which 
in turn is generated by the collective action of the atoms on the radiation (Prop. 4.1). 

(2) This microdynamics serves to drive the system into states in which the observables 
of different atoms are uncorrelated (Prop. 5.2). These states have particularly simple 
forms in both the normal and coherent phases, being stationary and space-translationally 
invariant in the former and periodic w.r.t. both space and time, in the latter. 

(3) In each phase, the eventual time-dependent microstate of the system maximises 
the specific entropy of the model, subject to the constraints imposed by the instantaneous 
values of the macroscopic variables. To the best of our knowledge, this is the first exam- 
ple where such a variational principle has been established in a regime far from thermal 
equilibrium. 

These results concern (1) an interplay between microphysics and macrophysics, (2) 
the decorrelation of certain micro-observables in the course of time, and (3) a maximum 
entropy principle, corresponding to a conditional thermodynamic stability far from thermal 
equilibrium. Although the model treated here is but a caricature of a laser, one might 
envisage that counterparts of these results might have at least some measure of validity 
in the physics of real dissipative systems. For example, one might expect a local version 
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of the maximum entropy principle (3) to prevail in hydrodynamics, in such a way as to 
impose local equilibrium conditions on the states of fluids. 
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